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Under the consideration of the solar energy potential of Turkey, a V-type Stirling engine 
having two heaters was designed, optimized and then manufactured. The prototype engine 
was tested in laboratory condition using an electrical heating system. Tests were con¬ 
ducted within the temperature range of 650-1000 °C with 50 °C increments. The pressure 
ranged from the ambient value to 2 bar with 0.5 bar increments at each stage of tempera¬ 
ture. The maximum power was obtained at 950 °C and 1.0 bar charge pressure as 118 W. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engines are externally heated, closed-cycle regenerative machines. Theoretical thermodynamic analysis of Stirling 
cycle engines indicates a thermal efficiency equal to that of the Carnot cycle. Although invented before internal combustion 
engines, Stirling engines are not yet commercially available today. The main reason preventing the Stirling engine from being 
commercial is its incompetitiveness with internal combustion engines in specific power and weight [1-3]. 

After inventing the Stirling engine in 1816, Robert Stirling carried out investigations up to 1850s and built double and 
triple cylinder engines, which were different in construction from the first engine. However, the performance of these en¬ 
gines was no better than that of the first engine. In 1860, Lehman built an horizontal single cylinder, displacer piston type 
Stirling engine without regenerator. This engine was found to be more popular than the regenerative engines built by Stirling 
[3,4]. In 1876, Rider used a novel construction on which expansion and compression spaces were set in different cylinders. 
On this engine, cylinders were present side by side and externally connected to each other by a regenerator. The superiority 
of this engine was the simple maintenance and silent operation as a result of having only five moving parts. Up to 1915, the 
only progress made on Stirling engines was the rearrangement of cylinders at right angles to each other by Robinson. The 
best Stirling engine manufactured by the end of 1923 had about 800 kg mass and 1.5 kW power output and 3% thermal effi¬ 
ciency [4,5]. 

In 1938, Professor Host noticed that there was a tremendous difference between the theoretical and the practical thermal 
efficiency of Stirling engines. Professor Host found this situation promising in consideration of the developments in material 
technology, and the research was started again in the Physical Research Laboratories in Eindhoven under the direction of H. 
Rinia. Several successful prototypes were built 50 times lighter in weight and 125 times smaller in volume than the engine 
designed in 1923. The results of their researches were published in 1947. These engines, called Philips Stirling Engines, 
worked under 20-50 bar charge pressure and had a maximum speed of 2000 rpm and an output of 2 HP [5,6]. 

In 1953, a new engine was invented by R.J. Meijer in Philips-Eindhoven, the Rhombic Drive Mechanism [7]. The primary 
advantage of this engine was that the piston applies no force to the cylinder wall and allows the use of a flexible seal between 
the piston and cylinder wall. Consequently, no charge pressure is needed in the crankcase. Meijer used hydrogen, helium and 
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air as working gases and obtained the best performance with hydrogen, the highest performance obtained up to that time 

[ 8 ]. 

After 1960, some engine companies (MAN-MWM, United Stirling of Sweden, Ford Motor Company of Detroit, etc.) started 
research programs to develop Stirling engines for automotive applications. Within these programs some multi-cylinder en¬ 
gines having different driving mechanisms were built and tested. The thermal efficiency of Stirling engines designed for 
automotive application was over 40%. However, the weights and volumes are not at desired levels [1]. 

In 1974, W. Beale invented the free-piston Stirling engine. In this engine, the phasing between working piston and dis¬ 
placer was not performed by a driving mechanism. Therefore, the engine had lower volume and weight, and thus, it had sev¬ 
eral advantages [9]. Particularly in solar energy research, it is appropriate to place the engine at the focus of a solar collector. 

The practical cycle of Stirling engines differs from the theoretical cycle consisting of two constant temperatures and two 
constant volume processes. For the thermodynamic analysis of Stirling engines, more authentic approximations were de¬ 
vised. Schmidt [10] introduced one of the approximations with the assumption that the fluids in hot and cold spaces are 
at the same temperature as the hot and cold sources, respectively. In this analysis instead of theoretical processes, the real 
processes performed by pistons moving sinusoidally with 90-degrees phase angle were used. Another thermodynamic anal¬ 
ysis was devised by Finkelstein [11], called nodal analysis, wherein the whole engine was divided into 13 sub-volumes. In 
this analysis, the temperatures at sub volumes are calculated with differential time intervals by the First Law of Thermody¬ 
namics. In the analysis, the instantaneous variations of the wall temperatures are also taken into account. Urieli et al. [12,13] 
improved the nodal analysis to cover the effects of leakage, heat losses and flow losses [14]. They conducted an investigation 
on nodal analysis, prepared computer programs and applied to the design of a rhombic driven engine for NASA. The details of 
nodal analysis can be found in their report. 

In the present study, a V-type Stirling engine prototype was designed and manufactured for solar energy applications and 
for domestic applications. The reserve of consumable energy sources such as oil and natural gas has been decreasing, and 
cost input for production has been increasing rapidly in recent days. It is expected that these sources will be run out within 
35-200 years [15]. In terms of solar energy, Turkey has a big potantice and solar energy is used for different purposes. In 
particular, during the summer time, solar energy per area is quite high [16]. 


2. Engine construction 

The engine was designed as V-(a) type with two heaters (double displaces). In V-type engine, there is a 90° phase angle 
between cylinders. The 90° phase angle enables the processes of the cycle. Schematic picture of Stirling engine is shown in 
Fig. 2.1. The area is shown in Fig. 2.1 as number 2 is called hot area for heating of process material and also number 6 is 
called cold area for cooling of process material. In order to alienate process material easily and to raise difference between 
cold and hot sources, cold area was cooled by cold water. The system is completely closed. Helium was used as a process 
material. The process material which is heated in number 2 hot area generates output by expanding in number 10 
cylinders. 

2 A. Heater and cooler 

Heater and cooler consist of four parts as displacer piston with center bearing, displacer piston, displacer cylinder, dis¬ 
placer connecting piston. Displacer piston provides to displace process materials (helium) by actuating between hot and cold 
areas, and meanwhile provide to transfer necessary heat energy to process materials with constant temperature in order to 
operate engine. Displacer piston operates in center bearing without lubrication. Center bearing was manufactured from hard 
cast iron with high graphite content in order to avoid inavodouble bearing friction. Displacer cylinder was manufactured 
from Cr alloyed steel with 3 mm thickness to rise its resistance against corrosion and heat. There is 1 mm gap between each 
displacer piston and cylinder. This gap also performs as regenerator. Heat is transferred from cylinder to process material by 
radiation. Whilst process material passes between piston and cylinder, its pressure and temperature is raised by absorbing 
heat and then output is generated. The temperature of displacer piston is less than that of cylinder surface and considering 
that while temperature of cylinder surface is at 800 K, temperature of displacer piston is 700 K. If the heat transferred by 
radiation is calculated between both surfaces, it is about 300 W. Some of this heat is transferred to process material and 
the rest is lost by different ways. 

2.2. Cylinders and piston 

Cylinder manufactured from high C cast iron is dry jacketed type. Cylinder jacket was driven into its hole with suitable 
tightness by pressing. 

Power piston generated output by consuming energy load in itself with process material. Meanwhile power piston is 
moved linearly. This action is transmitted to crankshaft by means of piston rod. Power piston was manufactured from high 
C cast iron without ring in order to reduce friction. An inertial force was very high due to piston weight. Piston weight was 
reduced to 300 g. Then engine speed was increased. A piston extension was attached to the piston top in order to minimize 
dead volumes in engine. 
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1- Diplacer Piston 

2- Heater 

3- Heat Pipe 

4- Inlet for cold water 

5- Outlet for cold water 

6- Cooling 

7- Centering bearing 

8- Connecting bar 

9- Crank shaft 

10- Power cylinder 

11- Power piston 

12- Power piston extension 

13- Power cylinder cover 

14- Block 



Fig. 2.1. Schematic drawing of manufactured Stirling engine with two displacer pistons. 


3. Thermodynamic analyses 

P-V and T-S diagrams of Stirling cycle are given in Fig. 3.1. Based on this cycle the operation of Stirling engines is ex¬ 
plained through Fig. 3.2 . The cycle consists of two isotherms (constant temperature) and two isochors (constant volume). 
The heat is given to process material which is placed in closed cylinder by means of heat exchanger heater. The exhausting 
of heat is performed by another heat exchanger cooler. Necessary heat for engine is provided from a special combustion 
chamber out of cylinder and a continuous combustion is maintained there. 



Fig. 3.1. P-V and 
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Fig. 3.2. Insiding (displacement) time. 
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Fig. 3.3. Stirling engine examined thermodynamically. 


In Stirling cycle first, heat is transferred to a process material along 3-4 isotherm at P 3 pressure and T H temperature then 
this provides for expansion up to P 4 pressure. Owing to taking heat from the system along 4-1 isochors latter on, the pres¬ 
sure of process material decreases to Pi and its temperature T c . After this operation, the extraction of heat from the system is 
maintained; process material is compressed along 1-2 isotherms up to P 2 . Finally, heat is given to system along 2-3 isochors, 
then, the temperature of process material is raised up to P 4 . 

3.1. Operation of Stirling engine 

Assuming that crankshaft moves from point A to point B, displacer piston is in mobility at this moment. Power piston 
moves from BDP to TDP and compresses process material in its front. This case is represented in PV diagram between num¬ 
bers 1 and 2. When crankshaft moves from B to C, displacer piston moves from top point to bottom point. It takes BC line in 
(At) time. This time is called siding (displacement). At this stage, it transfers process material from cold area to hot area. At 
this stage power piston is in mobility, waiting for change of direction. This is represented in PV diagram between 2 and 3. 

Process material with rising pressure and temperature provides crankshaft to move from C point to D point direction by 
impacting via power piston. At this stage output is obtained from the engine. Displacer piston is stable during this time and 
this is represented with points 3 and 4 in PV diagram. When crankshaft moves from point D to (A), displacer piston moves 
and it transfers process material from hot area to cold area. At this stage, the power piston is stable; this is represented be¬ 
tween points 4 and 1 in PV diagram. 

3.2. Analysis method 

The practical Stirling cycle consists of two isothermal and two isochor processes, rather different from ideal Stirling cycle. 
For this reason, it is used as rather realistic approach to calculate thermodynamic values during designing stage. 

In this study design details of the manufactured engine were obtained by using simulation software developed by Kar- 
abulut. Thermodynamic simulation software is given in Annex. F. the designed engine has five volumes and thermodynam¬ 
ically investigated areas are schematically shown in Figs. 3-4. 

Following assumptions are used in analysis method: 

1. Every cell is an open system from the thermodynamic point of view and the cell has material inlet and outlet under the 
periodic conditions. 

2. Pressure differences which arise from hydrodynamic frictions during the displacement of process material can be 
neglected comparing to thermodynamic pressure. 
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3. There is a linear and stable temperature distribution in the stiff part of regenerator and it provides heat balance for the 
stiff part of regenerator. 

4. Helium is considered ideal gas and used as a process material. 

5. Total mass in the system does not change with time. 

6. The walls of hot cylinder and heater are at hot source temperature and the walls of cold cylinder and cooler are at cold 
source temperature. 

7. The volumes in heater, cooler and regenerator are constant; the volumes in hot and cold cylinders alter according to the 
crank angle, and they can be expressed by following formula: 


V c = V D +^Vs[l+cos(8-<!>)] 
Vh = V D + ^ V s [l + cos(0)] 


(1) 

( 2 ) 


There are three unknown terms such as pressure, temperature and density and the formula which can be used to deter¬ 
mine the above terms can be given as follows based on formula for ideal gases, mass balance and energy balance, 
respectively: 

P = PRT (3) 

m A + m 2 + m 3 h -b m n = m t (4) 

dQ - dW + dHjniet - dHoutiet = dU (5) 


Taking into consideration that pressure is the same at every condition and (3) and (4) can be evaluated as follows: 


V 


m t R 


Tl I ^2 I I | V n 

T-[ ' t ' t 3 ' h r„ 



The formula given below can be arranged according to energy balance: 

hiAi(T wi - Ti)dt - pdVi + (dH in i et - dH 0U tiet) = lriiCydTi + C v Tidmi (7) 

where (dH in i et - dH ou tiet) term in the formula presents input and output energy from component bar dens and is calculated as 
follows: 


(dH in iet - dHoutiet)* = C p [(m\ + + • • • + m\_ A ) - (m\ +At + rri+ At + 


+ mf+f)] 


Ti -1 
Ti 


+ Cp[( m i+1 + m i+2 ^-f m n) 


- (mj+f + At -i-h m^ +Af )] 


Ti 

Tt+i 


In the above formula, the first term shows the energy which enters or leave from the component from the left side; the 
second term shows the energy which enter to or leave from the component from the right side. The temperatures are written 
in matrix column; those are located upper are used for flow which moves from left to right and those located lower are used 
for flow which moves from right to left. The formula used numerical procedures as follows [17]: 

AT, = - T,)At - (m‘ +At - m‘)C v m‘)C v T i+ 

(dHiniet - dH 0Utl et)i - PAVi]/C v m, 


P-V diagram which is designed and obtained by software by using this analysis method is shown in Fig. 3.4. 



Fig. 3.4. P-V diagram obtained by software. 
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Table 3.1 

Calculated and measured engines values 


Properties 

Obtained by software 


Type of engine 

V (a) 


Hot cylinder volume (V E ) 

311 cm 3 


Cold cylinder volume (V H ) 

311 cm 3 


Total dead volume (V D ) 

290.57 cm 3 


Diameter of cylinder (D) 

61 mm 


Piston stroke (L) 

65 mm 


Fazes angle (</>) 

90° 


Cold source temperature (T c ) 

350 I< 

318 K 

Hot source temperature (T H ) 

1273 K 

953 K 

Maximum engine speed (n max ) 

1200 1/min 

784 1/min 

Carnot cycle (tjQ 

% 72 a 

% 0.66 

Net work (W net ) 

245 J 


Engine efficiency (//) 

11% 


Charge pressure (P C harge), (manometric) 

2 bar 

2 bar 

Max cycle pressure (P max ) 

3.048 bar 3 


Mean cycle pressure (P ort ) 

2.810 bar 3 


Process material 

Helium 


Output (N) 

500 W 

118 W 

Compression ratio (e) 

1.65 3 


Minimum volume (V min ) 

290.57 cm 3 


Maximum volume (V max ) 

480.53 cm 3 



a Value based on theoretical calculation. 
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Fig. 3.5. Variation between gas and wall temperatures depending on the length of displacer piston. 


The values belonging to Stirling engine by using software are shown is Table 3.1. 

The manufactured engine was designed for 500 W output. However 118 W of output was obtained during the tests. One 
of the reasons for low output is access section distance (( 5 ). Since the access section gets more narrow and since gas temper¬ 
ature gets closer to wall temperature, engine efficiency will rise, if the distance is far away as much as possible between hot 
and cold source according to Carnot cycle. On the contrary, the more the access section distance rises, the less the differences 
between cold and hot sources, hence the efficiency gets lower. 

Different access section with various gas temperatures and wall temperatures are shown in Fig. 3.5 based on the length of 
the displacer piston. As it is seen in this figure, wall temperature and gas temperature vary linearly depending on the length 
of the displacer piston. However, since the distance between displacer piston and cylinder (access sections) reduces, the bet¬ 
ter heat transmission is obtained; gas temperature gets closer to wall temperature. On the other hand, even it is clear that 
when the access section decreases more, consequently heat transmission gets well, since access of process material becomes 
difficult; losing output (pumping work) gets increased. When taking manufacturing into account, it is difficult to centralize 
displacer piston, if the access section is very small. If the access section distance is very high, heat transmission gets worse; 
therefore, gas temperature recedes from wall temperature. The access section is given 1 mm width in manufacturing. 

4. Engine performance 

Test engine was operated at constant oven temperatures without charge, and the change of speed, torque, output vibra¬ 
tion was measured in the same conditions, but at the different charge pressure speed, torque, output and vibration were 
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Fig. 4.1. Variation of engine speed depending on heater temperature (without charge). 


examined (gauge pressures of 0.5 bar, 1.0 bar, 1.5 bar, 2.0 bar). Fig. 4.1 shows the effect of heater temperature on engine 
speed. 

Variation of engine speed depending on heater temperature is shown in Fig. 4.1. As it is seen in this figure, while heater 
temperature increases, engine speed also increases linearly. Engine speed varies based on heater temperature. In Fig. 4.2 
variations of engine speed depending on heater temperature are shown at different charge pressures. As it is seen in this 
figure, the engine speed of 784 1/min was reached at 10 bar pressure, 1000 °C, and speeds vary depending on heater temper¬ 
ature. Engine speed of 645 1/min was reached under the atmospheric condition at 1000 °C (Fig. 4.1). The engine speed was 
high at charge conditions compared to atmospheric conditions, owing to the increase of mass of process materials. 

4.1. Impact of heater temperature on torque and output 

The variation of engine torque and output depending on heater temperature is shown in Fig. 4.3, output of 34 W was ob¬ 
tained at 950 °C. A graph was created by taking maximum torque and output values for each heater temperature; increasing 
of torque and output depending on heater temperatures causes different temperatures between heater and cooler. Owing to 
heater temperature increase up to material endurance limit raises cycle output and engine efficiency, torque and output also 
rise depending on it. 

The variation of engine torque and output at different charge pressures and heater temperatures is shown Fig. 4.4. Max¬ 
imum power of 118 W was obtained at 950 °C and at 1.0 bar pressure. When looked at torque curve, the same situation as 
power curve can be seen . The reason for decrease of torque and output after 1.0 bar is due to the inadequacy of temperature 
heat transfer surface area. 

5. Discussion and outlook 

Engine efficiency was 11%. It is necessary to eliminate dead volumes and leakage in the system which increases engine 
efficiency. 

Engine speed increases depending on heater temperature. Increasing of heater temperature up to material endurance lim¬ 
its was very effective in improving engine torque and output. 



Fig. 4.2. Variation of engine speed depending on heater temperature (with charge). 
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600 700 800 900 1000 


Heater Temperature (C) 

Fig. 4.3. Variation of torque and output depending on heater temperature (without charge). 



600 700 800 900 1 000 

Heater Temperature (C) 

Fig. 4.4. Variation of torque and output depending on heater temperature (at 1 bar pressure). 

Charge pressure to system was performed from carter. In order to decrease the leakage, it will be high efficiency, if the 
charge performs only from the heater and chiller sections. 

In this study, it was designed to obtain 500 W of output, but output of 118 W was obtained at 594 1/min during the tests. 
In order to take more output, it is necessary to minimize dead volumes, to well balance between crank and rod mechanism, 
and provide very good leakage to charge the system. 

Another use which needs to solve is to make Stirling engine compete with internal combustion engine rather than leak¬ 
age, dead volumes and improvement of heat transmission, storage of solar energy. Due to solar energy being familiar to envi¬ 
ronment and very cheaper, same collectors which are used to collect solar energy can be used so that Stirling engines can be 
used for irrigation purpose. 

Manufactured V-type engine with double displacer pistons is an example which is known in Turkey. Prototype balancing, 
leakage, heat transmission and dead volumes still exist. If these problems are solved, they can be used as stationary plant 
engine in the market. 
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